A framework based on subject-specific musculoskeletal models and Monte Carlo simulations to personalize muscle coordination retraining.

Excessive loads at lower limb joints can lead to pain and degenerative diseases. Altering joint loads with muscle coordination retraining might help to treat or prevent clinical symptoms in a non-invasive way. Knowing how much muscle coordination retraining can reduce joint loads and which muscles have the biggest impact on joint loads is crucial for personalized gait retraining. We introduced a simulation framework to quantify the potential of muscle coordination retraining to reduce joint loads for an individuum. Furthermore, the proposed framework enables to pinpoint muscles, which alterations have the highest likelihood to reduce joint loads. Simulations were performed based on three-dimensional motion capture data of five healthy adolescents (femoral torsion 10°-29°, tibial torsion 19°-38°) and five patients with idiopathic torsional deformities at the femur and/or tibia (femoral torsion 18°-52°, tibial torsion 3°-50°). For each participant, a musculoskeletal model was modified to match the femoral and tibial geometry obtained from magnetic resonance images. Each participant's model and the corresponding motion capture data were used as input for a Monte Carlo analysis to investigate how different muscle coordination strategies influence joint loads. OpenSim was used to run 10,000 simulations for each participant. Root-mean-square of muscle forces and peak joint contact forces were compared between simulations. Depending on the participant, altering muscle coordination led to a maximum reduction in hip, knee, patellofemoral and ankle joint loads between 5 and 18%, 4% and 45%, 16% and 36%, and 2% and 6%, respectively. In some but not all participants reducing joint loads at one joint increased joint loads at other joints. The required alteration in muscle forces to achieve a reduction in joint loads showed a large variability between participants. The potential of muscle coordination retraining to reduce joint loads depends on the person's musculoskeletal geometry and gait pattern and therefore showed a large variability between participants, which highlights the usefulness and importance of the proposed framework to personalize gait retraining.

Comprehensively characterizing heterogeneous and transversely isotropic properties of femur cortical bones.

Comprehensive characterization of the transversely isotropic mechanical properties of long bones along both the longitudinal and circumferential gradients is crucial for developing accurate mathematical models and studying bone biomechanics. In addition, mechanical testing to derive elastic, plastic, and failure properties of bones is essential for modeling plastic deformation and failure of bones. To achieve these, we machined a total of 336 cortical specimens, including 168 transverse and 168 longitudinal specimens, from four different quadrants of seven different sections of 3 bovine femurs. We conducted three-point bending tests of these specimens at a loading rate of 0.02 mm/s. Young's modulus, yield stress, tangential modulus, and effective plastic strain for each specimen were derived from correction equations based on classical beam theory. Our statistical analysis reveals that the longitudinal gradient has a significant effect on the Young's modulus, yield stress, and tangential modulus of both longitudinal and transverse specimens, whereas the circumferential gradient significantly influences the Young's modulus, yield stress, and tangential modulus of transverse specimens only. The differences in Young's modulus and yield stress between longitudinal specimens from different sections are greater than 40%, whereas those between transverse specimens are approximately 30%. The Young's modulus and yield stress of transverse specimens in the anterior quadrant were 18.81%/15.46% and 18.34%/14.88% higher than those in the posterior and lateral quadrants, respectively. There is no significant interaction between the longitudinal gradient and the circumferential gradient. Considering the transverse isotropy, it is crucial to consider loading direction when investigating the impact of circumferential gradients in the anterior, lateral, medial, and posterior directions. Our findings indicate that the conventional assumption of homogeneity in simulating the cortical bone of long bones may have limitations, and researchers should consider the anatomical position and loading direction of femur specimens for precise prediction of mechanical responses.

Muscle forces acting on the greater trochanter lead to a dorsal warping of the apophyseal growth plate.

The apophyseal growth plate of the greater trochanter, unlike most other growth plates of the human body, exhibits a curved morphology that results in a divergent pattern resembling an open crocodile mouth on plain antero-posterior radiographs. To quantify the angular alignment of the growth plate and to draw conclusions about the function of the muscles surrounding it, we analyzed 57 MRI images of 51 children and adolescents aged 3-17 years and of six adults aged 18-52 years. We measured the angulation of the plate relative to the horizontal plane (AY angle) and the trajectories of the muscles attaching to the greater trochanter of the proximal femur. From anterior to posterior, the AY angle shows a decrease of 33.44°. In the anterior third, the cartilage is angled at a mean of 51.64°, and in the posterior third, the mean angulation is 18.6°. This indicates that the cartilage in the anterior region of the greater trochanteric apophysis is subject to more vertically oriented force vectors compared to the posterior region, as the growth plates align perpendicular to the force vectors acting on them. Combining the measured muscle trajectories with the physiological cross-sectional areas (PCSA) available from the literature revealed that, in addition to the known internal and external lateral traction ligament systems, a third, dorsally located traction ligament system exists that may be responsible for the dorsal deformation of the AY angle.

Influence of femoral anteversion angle and neck-shaft angle on muscle forces and joint loading during walking.

Femoral deformities, e.g. increased or decreased femoral anteversion (AVA) and neck-shaft angle (NSA), can lead to pathological gait patterns, altered joint loads, and degenerative joint diseases. The mechanism how femoral geometry influences muscle forces and joint load during walking is still not fully understood. The objective of our study was to investigate the influence of femoral AVA and NSA on muscle forces and joint loads during walking. We conducted a comprehensive musculoskeletal modelling study based on three-dimensional motion capture data of a healthy person with a typical gait pattern. We created 25 musculoskeletal models with a variety of NSA (93°-153°) and AVA (-12°-48°). For each model we calculated moment arms, muscle forces, muscle moments, co-contraction indices and joint loads using OpenSim. Multiple regression analyses were used to predict muscle activations, muscle moments, co-contraction indices, and joint contact forces based on the femoral geometry. We found a significant increase in co-contraction of hip and knee joint spanning muscles in models with increasing AVA and NSA, which led to a substantial increase in hip and knee joint contact forces. Decreased AVA and NSA had a minor impact on muscle and joint contact forces. Large AVA lead to increases in both knee and hip contact forces. Large NSA (153°) combined with large AVA (48°) led to increases in hip joint contact forces by five times body weight. Low NSA (108° and 93°) combined with large AVA (48°) led to two-fold increases in the second peak of the knee contact forces. Increased joint contact forces in models with increased AVA and NSA were linked to changes in hip muscle moment arms and compensatory increases in hip and knee muscle forces. Knowing the influence of femoral geometry on muscle forces and joint loads can help clinicians to improve treatment strategies in patients with femoral deformities.

The application of an isotropic crushable foam model to predict the femoral fracture risk.

For biomechanical simulations of orthopaedic interventions, it is imperative to implement a material model that can realistically reproduce the nonlinear behavior of the bone structure. However, a proper material model that adequately combines the trabecular and cortical bone response is not yet widely identified. The current paper aims to investigate the possibility of using an isotropic crushable foam (ICF) model dependent on local bone mineral density (BMD) for simulating the femoral fracture risk. The elastoplastic properties of fifty-nine human femoral trabecular cadaveric bone samples were determined and combined with existing cortical bone properties to characterize two forms of the ICF model, a continuous and discontinuous model. Subsequently, the appropriateness of this combined material model was evaluated by simulating femoral fracture experiments, and a comparison with earlier published results of a softening Von-Mises (sVM) material model was made. The obtained mechanical properties of the trabecular bone specimens were comparable to previous findings. Furthermore, the ultimate failure load predicted by the simulations of femoral fractures was on average 79% and 90% for the continuous and discontinuous forms of the ICF model and 82% of the experimental value for the sVM material model. Also, the fracture locations predicted by ICF models were comparable to the experiments. In conclusion, a nonlinear material model dependent on BMD was characterized for human femoral bone. Our findings indicate that the ICF model could predict the femoral bone strength and reproduce the variable fracture locations in the experiments.

Investigation of Characteristic Motion Patterns of the Knee Joint During a Weightbearing Flexion.

This study aimed to develop and validate a novel flexion axis concept by calculating the points on femoral condyles that could maintain constant heights during knee flexion. Twenty-two knees of 22 healthy subjects were investigated when performing a weightbearing single leg lunge. The knee positions were captured using a validated dual fluoroscopic image system. The points on sagittal planes of the femoral condyles that had minimal changes in heights from the tibial plane along the flexion path were calculated. It was found that the points do formulate a medial-lateral flexion axis that was defined as the iso-height axis (IHA). The six degrees of freedom (6DOF) kinematics data calculated using the IHA were compared with those calculated using the conventional transepicondylar axis and geometrical center axis. The IHA measured minimal changes in proximal-distal translations and varus-valgus rotations along the flexion path, indicating that the IHA may have interesting clinical implications. Therefore, identifying the IHA could provide an alternative physiological reference for improvement of contemporary knee surgeries, such as ligament reconstruction and knee replacement surgeries that are aimed to reproduce normal knee kinematics and medial/lateral soft tissue tensions during knee flexion.

Stromal Lineage Precursors from Rodent Femur and Tibia Bone Marrows after Hindlimb Unloading: Functional Ex Vivo Analysis.

Rodent hindlimb unloading (HU) model was developed to elucidate responses/mechanisms of adverse consequences of space weightlessness. Multipotent mesenchymal stromal cells (MMSCs) were isolated from rat femur and tibia bone marrows and examined ex vivo after 2 weeks of HU and subsequent 2 weeks of restoration of load (HU + RL). In both bones, decrease of fibroblast colony forming units (CFU-f) after HU with restoration after HU + RL detected. In CFU-f and MMSCs, levels of spontaneous/induced osteocommitment were similar. MMSCs from tibia initially had greater spontaneous mineralization of extracellular matrix but were less sensitive to osteoinduction. There was no recovery of initial levels of mineralization in MMSCs from both bones during HU + RL. After HU, most bone-related genes were downregulated in tibia or femur MMSCs. After HU + RL, the initial level of transcription was restored in femur, while downregulation persisted in tibia MMSCs. Therefore, HU provoked a decrease of osteogenic activity of BM stromal precursors at transcriptomic and functional levels. Despite unidirectionality of changes, the negative effects of HU were more pronounced in stromal precursors from distal limb-tibia. These observations appear to be on demand for elucidation of mechanisms of skeletal disorders in astronauts in prospect of long-term space missions.

The effects of orthodontic anchor screw inserted into the femur of growth-phase or mature rats -Osteoid formation, bone mineral density, collagen fiber bundles, biological apatite crystal orientation.

The purpose of this study is to investigate the effect of orthodontic anchor screws (OASs) inserted into the femur of growth-phase or mature rats using histological observation and bone structure analysis. The experimental animals are growth-phase (6-week-old) or mature (25-week-old) male Wistar rats. OAS was placed into the point one-third of the femoral length from the proximal end of the femur, and the response of the surrounding bone was observed and measured. The results showed at the OAS bone interface, in growth-phase rats, bone mineral density (BMD) was reduced and the running angle of collagen fiber bundles varied significantly. In mature rats, more osteoid was observed and biological apatite (BAp) crystals showed a different orientation. It was suggested that after the insertion of OASs, bone volume and quality are decreased, but after a sufficient healing period, a new bone micro/nano structure, different from the original structure, are reconstructed.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy.

BACKGROUND AND OBJECTIVE: Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure. METHODS: hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases. RESULTS: The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role. CONCLUSIONS: hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.

Normal coronal kinematics of dynamic alignment and bony positions relative to the ground in three-dimensional motion analysis during gait: A preliminary study.

BACKGROUND: During gait, healthy knee coronal kinematics of each bony axis and lower extremity alignment are important because they could be useful as reference data for several surgeries and provide clarification of the etiology of diseases around the knee in healthy participants; however, it remains unknown. OBJECTIVE: The objective of this study was to clarify the kinematics of lower extremity alignment and the bony axes relative to the ground during gait, focused on the coronal plane, in healthy individuals by applying our unique three-dimensional (3D) motion analysis. METHODS: The study included 21 healthy individuals, including 9 healthy females and 12 healthy males with an average age of 36 ± 17 years. Knee kinematics were calculated in a gait analysis by combining the data from a motion-capture system and a 3D lower-extremity alignment assessment system on biplanar long-leg radiographs by using a 3D-2D registration technique. The main kinematic parameters were the dynamic position change relative to the ground, applying the femoral anatomical axis (FAA), tibial anatomical axis (TAA), and dynamic alignment in the coronal plane during the stance phase of gait. RESULTS: The average changes in FAA, TAA, and dynamic varus alignment were 3.7° ± 1.2°, 3.5° ± 0.8°, and 3.0° ± 1.2°, respectively. The TAA tilted laterally during the loading response and a plateau area appeared afterwards; the FAA gradually inclined laterally until the terminal stance phase, and the dynamic alignment showed varus angular change during the loading response. CONCLUSIONS: The tibia and femur were found to change approximately 2-5° of the position of the bony axes relative to the ground. In terms of clinical relevance, our findings can be used to clarify the etiology of diseases around the knee joint and as reference data for surgeries.

Computational prediction of the long-term behavior of the femoral density after THR using the Silent Hip stem.

Aseptic loosening due to the progressive periprosthetic bone resorption following total hip replacement is a crucial concern, that causes complications and failure of the arthroplasty surgery. The mismatch in stiffness between the hip implant and the surrounding femoral bone is one of the key factors leading to bone density resorption. This paper aimed to investigate the long-term response of the femoral bone after THR using the Silent Hip stem. For this purpose, a validated thermodynamic-based computational model was used to compute the change in bone density before and after THR. This model incorporated essential factors involved in bone remodeling process, such as mechanical loading, and biochemical affinities. The results of the numerical simulations using 3D finite element analysis were analyzed in five zones of interest qualitatively and quantitatively. Bone density predictions showed notable bone resorption in cervical areas, specifically in zone 1 and zone 5 of -18.7% and -14%, respectively. Conversely, bone formation was observed in the greater trochanter area (zone 2) of +25%. Stress shielding seemed to occur at cervical area due to the reduction in the mechanical loading in this region. Based on the quantitative analysis of the bone density distribution throughout the femoral bone, it appears that the Silent Hip stem achieved less bone resorption compared to conventional hip stem designs reported in the literature, which could be used for active patients.

A unified framework of cell population dynamics and mechanical stimulus using a discrete approach in bone remodelling.

Multiphysics models have become a key tool in understanding the way different phenomenon are related in bone remodeling and various approaches have been proposed, yet, to the best of the author's knowledge there is no model able to link a cell population model with a mechanical stimulus model using a discrete approach, which allows for an easy implementation. This article couples two classical models, the cell population model from Komarova and the Nackenhorst model in a 2D domain, where correlations between the mechanical loading and the cell population dynamics can be established, furthermore the effect of different paracrine and autocrine regulators is seen on the overall density of a portion of trabecular bone. A discretization is performed using frame 1D finite elements, representing the trabecular structure. The Nackenhorst model is implemented by using the finite element method to calculate the strain energy as the main mechanical stimulus that determines the bone mass density evolution in time. This density is normalized to be added to the bone mass percentage proposed by the Komarova model, where coupling terms have been added as well that guarantee a stable response. In the simulations, the equations were solved employing the finite element method with a user subroutine implemented in ABAQUS (2017) and by applying a direct formulation. The methodology presented can model the cell dynamics occurring in bone remodelling in accordance with the asynchronous nature of this process, yet allowing to differentiate zones with higher density, the main trabecular groups are obtained for the proximal femur. Finally, the model is tested in pathological cases, such as osteoporosis and osteopetrosis, yielding results similar to the pathology behavior. Furthermore, the discrete modelling technique is shown to be of use in this particular application.

Physiological Axial Tibial Rotation of the Knee During a Weightbearing Flexion.

Axial tibial rotation is a characteristic motion of the knee, but how it occurs with knee flexion is controversial. We investigated the mechanisms of tibial rotations by analyzing in vivo tibiofemoral articulations. Twenty knees of 20 living human subjects were investigated during a weightbearing flexion from full extension to maximal flexion using a dual fluoroscopic imaging system. Tibiofemoral articular contact motions on medial and lateral femoral condyles and tibial surfaces were measured at flexion intervals of 15 deg from 0 deg to 120 deg. Axial tibial rotations due to the femoral and tibial articular motions were compared. Articular contact distances were longer on femoral condyles than on tibial surfaces at all flexion intervals (p < 0.05). The articular distance on medial femoral condyle is longer than on lateral side during flexion up to 60 deg. The internal tibial rotation was 6.8 ± 4.5 deg (Mean ± SD) at the flexion interval of 0-15 deg, where 6.1 ± 2.6 deg was due to articulations on femoral condyles and 0.7 ± 5.1 deg due to articulations on tibial surfaces (p < 0.05). The axial tibial rotations due to articulations on femoral condyles are significantly larger than those on tibial surfaces until 60 deg of flexion (p < 0.05). Minimal additional axial tibial rotations were observed beyond 60 deg of flexion. The axial tibial rotations were mainly attributed to uneven articulations on medial and lateral femoral condyles. These data can provide new insights into the understanding of mechanisms of axial tibial rotations and serve as baseline knowledge for improvement of knee surgeries.

Effects of hip muscle activation on the stiffness and energy absorption of the trochanteric soft tissue during impact in sideways falls.

The trochanteric soft tissue attenuates impact force or absorbs impact energy during a fall on the hip (thereby helps to reduce a risk of hip fracture). While the benefits should be affected by contractions of muscles spanning the hip joint, no information is available to date. We examined how the stiffness (force attenuation capacity) and energy absorption of the trochanteric soft tissue were affected by hip muscle activation during a fall. Thirteen healthy young individuals (5 males, 8 females) participated in the pelvis release experiment. Falling trials were acquired with three muscle contraction conditions: 0-20% ("relaxed"), 20-50% ("moderate"), and 60-100% ("maximal") of the maximal voluntary isometric contraction of the gluteus medius muscle. During trials, we measured real-time force and deformation behaviour of the trochanteric soft tissue. Outcome variables included the stiffness and energy absorption of the soft tissue. The stiffness and energy absorption ranged from 56.1 to 446.9 kN/m, and from 0.15 to 2.26 J, respectively. The stiffness value increased with muscle contraction, and 59% greater in "maximal" than "relaxed" condition (232.2 (SD = 121.4) versus 146.1 (SD = 49.9)). However, energy absorption decreased with muscle contraction, and 58.9% greater in "relaxed" than "maximal" condition (0.89 (SD = 0.63) versus 0.56 (SD = 0.41)). Our results provide insights on biomechanics of the trochanteric soft tissue ("natural" padding device) during impact stage of a fall, suggesting that soft tissues' protective benefits are largely affected by the level of muscle contraction.

Nordic walking with an integrated resistance shock absorber affects the femur strength and muscles torques in postmenopausal women.

Deterioration of the structure and function of the musculoskeletal system represents a significant problem during aging and intervention with a suitable load of physical activity may improve the quality of life. Nordic walking (NW) has become a popular and easily accessible form of activity, especially for older adults people around the world. Thus, the purpose of the study was to evaluate the influence of an Nordic walking training program with classic poles (NW) and with integrated resistance shock absorber (RSA) on bone mineral density and the peak torques of upper limb muscles and to compare the effects of both intervention programs. 25 women were randomly assigned to two training groups: 10 subjects using RSA (68 ± 4.19 years) and 15 subjects using NW poles (65 ± 3.40 years), which completed 8 weeks of training program. The hip, spine and forearm areal bone mineral density, torques of the flexors and extensors at the elbow and shoulder joints were measured before starting the training programs and after their completion. The most significant effect was found in differences between the two groups of women with respect to the femur strength index (p = 0.047) and the ratio of the flexors to extensors in the elbow (p = 0.049) and shoulder (p = 0.001) joints and peak torque of flexors in the shoulder joint (p = 0.001) for the left arm. A significant difference was also found in the index of torque asymmetry of flexors in the shoulder joint (p = 0.002). The study shows that Nordic walking with RSA poles for postmenopausal women led to beneficial changes in the femur strength index. However, we found no significant influence on bone mineral density values measured on the whole body, the femoral neck, forearm or lumbar spine regions. The occurrence of asymmetry in biomechanical muscle parameters, which was observed using RSA poles, may suggest the necessity of systematic controlling the gait technique to avoid the adverse consequences of asymmetrical rotation of the lumbar spine.

Dietary Hempseed Decreases Femur Maximum Load in a Young Female C57BL/6 Mouse Model but Does Not Influence Bone Mineral Density or Micro-Architecture.

Numerous seed and seed extract diets have been investigated as a means of combating age-related bone loss, with many findings suggesting that the seeds/extracts confer positive effects on bone. Recently, there has been rising interest in the use of dietary hempseed in human and animal diets due to a perceived health benefit from the seed. Despite this, there has been a lack of research investigating the physiologic effects of dietary hempseed on bone. Previous studies have suggested that hempseed may enhance bone strength. However, a complete understanding of the effects of hempseed on bone mineralization, bone micro-architecture, and bone biomechanical properties is lacking. Using a young and developing female C57BL/6 mouse model, we aimed to fill these gaps in knowledge. From five to twenty-nine weeks of age, the mice were raised on either a control (0%), 50 g/kg (5%), or 150 g/kg (15%) hempseed diet (n = 8 per group). It was found that the diet did not influence the bone mineral density or micro-architecture of either the right femur or L5 vertebrae. Furthermore, it did not influence the stiffness, yield load, post-yield displacement, or work-to-fracture of the right femur. Interestingly, it reduced the maximum load of the right femur in the 15% hempseed group compared to the control group. This finding suggests that a hempseed-enriched diet provides no benefit to bone in young, developing C57BL/6 mice and may even reduce bone strength.

Hind limb muscles influence the architectural properties of long bones in frogs.

The Mechanostat Theory states that osteocytes sense both the intensity and directionality of the strains induced by mechanical usage and modulate the bone design accordingly. In long bones, this process may adapt anterior-posterior and lateral-medial strength to their mechanical environment showing regional specificity. Anuran species are ideal for analyzing the muscle-bone relationships related to the different mechanical stresses induced by their many locomotor modes and habitat uses. This work aimed to explore the relationships between indicators of the force of the most relevant muscles to locomotion and the mechanical properties of femur and tibia fibula in preserved samples of three anuran species with different habitat use (aquatic, arboreal) and locomotion modes (swimmer, jumper, walker/climber). For that purpose, we measured the anatomical cross-sectional area of each dissected muscle and correlated it with the moments of inertia and bone strength indices. Significant, species-specific covariations between muscle and bone parameters were observed. Pseudis platensis, the aquatic swimmer, showed the largest muscles, followed by Boana faber, the jumper and Phyllomedusa sauvagii, the walker/climber. As we expected, bigger muscles correlate with bone parameters in all the species. Nevertheless, smaller muscles also play an important role in bone design. In aquatic species, muscle interaction enhances mostly lateral bending strength throughout the femur and lateral and antero-posterior bending strength in the tibia fibula. In the jumper species, muscles affected the femur and tibia fibula mostly in anterior-posterior bending. In the walker/climber species, responses involving both antero-posterior and lateral bending strengths were observed in the femur and tibia fibula. These results show that bones will be more or less resistant to lateral and antero-posterior bending according to the different mechanical challenges of locomotion in aquatic vs. arboreal habitats. This study provides new evidence of the muscle-bone relationships in three frog species associated with their different locomotion and habitat uses, highlighting the crucial role of muscle in determining the architectural properties of bones.

Determination of muscle strength and function in plesiosaur limbs: finite element structural analyses of Cryptoclidus eurymerus humerus and femur.

Background: The Plesiosauria (Sauropterygia) are secondary marine diapsids. They are the only tetrapods to have evolved hydrofoil fore- and hindflippers. Once this specialization of locomotion had evolved, it remained essentially unchanged for 135 Ma. It is still controversial whether plesiosaurs flew underwater, rowed, or used a mixture of the two modes of locomotion. The long bones of Tetrapoda are functionally loaded by torsion, bending, compression, and tension during locomotion. Superposition of load cases shows that the bones are loaded mainly by compressive stresses. Therefore, it is possible to use finite element structure analysis (FESA) as a test environment for loading hypotheses. These include muscle reconstructions and muscle lines of action (LOA) when the goal is to obtain a homogeneous compressive stress distribution and to minimize bending in the model. Myological reconstruction revealed a muscle-powered flipper twisting mechanism. The flippers of plesiosaurs were twisted along the flipper length axis by extensors and flexors that originated from the humerus and femur as well as further distal locations. Methods: To investigate locomotion in plesiosaurs, the humerus and femur of a mounted skeleton of Cryptoclidus eurymerus (Middle Jurassic Oxford Clay Formation from Britain) were analyzed using FE methods based on the concept of optimization of loading by compression. After limb muscle reconstructions including the flipper twisting muscles, LOA were derived for all humerus and femur muscles of Cryptoclidus by stretching cords along casts of the fore- and hindflippers of the mounted skeleton. LOA and muscle attachments were added to meshed volumetric models of the humerus and femur derived from micro-CT scans. Muscle forces were approximated by stochastic iteration and the compressive stress distribution for the two load cases, "downstroke" and "upstroke", for each bone were calculated by aiming at a homogeneous compressive stress distribution. Results: Humeral and femoral depressors and retractors, which drive underwater flight rather than rowing, were found to exert higher muscle forces than the elevators and protractors. Furthermore, extensors and flexors exert high muscle forces compared to Cheloniidae. This confirms a convergently evolved myological mechanism of flipper twisting in plesiosaurs and complements hydrodynamic studies that showed flipper twisting is critical for efficient plesiosaur underwater flight.

Stochastic failure analysis of proximal femur using an isogeometric analysis based nonlocal gradient-enhanced damage model.

BACKGROUND AND OBJECTIVE: Medical imaging-based finite element methods are more accurate tools for fracture risk prediction than the traditional aBMD based methods. However, these methods have drawbacks like geometric errors, high computational cost, mesh-dependent results, etc. In this article, the authors have proposed an isogeometric analysis-based nonlocal gradient-enhanced damage model to overcome some of these issues. Moreover, there are uncertainties in the values of input parameters for such analysis due to various measurement errors. Hence, stochastic analysis is performed to quantify the effect of these parametric uncertainties on the fracture behavior of the proximal femur. METHODS: Computed Tomography images of a patient are used to create a 2D proximal femur model with a heterogeneous description of material properties. A numerical model based on gradient-enhanced nonlocal continuum damage mechanics is used for fracture analysis of proximal femur to overcome the issues related to mesh dependency in traditional continuum damage mechanics models. Further, a multipatch isogeometric solver is developed to solve the governing equations. Monte Carlo simulations are used to understand the effect of parametric uncertainties on the fracture behavior of the proximal femur. RESULTS: The developed numerical framework is used to solve the fracture problem of proximal femur under single leg stance loading conditions. The obtained results are validated by comparing the load-displacement response and the crack path with that given in the literature. Stochastic analysis is performed by considering a ±5% variation in the elastic modulus, damage initiation strain, and the neck-shaft angle values. CONCLUSION: The proposed numerical framework can correctly predict the damage initiation and propagation in a proximal femur. The results reveal that the heterogeneous nature of material properties of bone plays a significant role in determining the fracture characteristics of the proximal femur. Further, the results of the stochastic analysis reveal that the parametric uncertainties in the neck-shaft angle have a much more significant influence on the results of the analysis than the parametric uncertainties in the elastic modulus and damage initiation strain.

ISB clinical biomechanics award winner 2021: Tibio-femoral kinematics of natural versus replaced knees - A comparison using dynamic videofluoroscopy.

BACKGROUND: A comparison of natural versus replaced tibio-femoral kinematics in vivo during challenging activities of daily living can help provide a detailed understanding of the mechanisms leading to unsatisfactory results and lay the foundations for personalised implant selection and surgical implantation, but also enhance further development of implant designs towards restoring physiological knee function. The aim of this study was to directly compare in vivo tibio-femoral kinematics in natural versus replaced knees throughout complete cycles of different gait activities using dynamic videofluoroscopy. METHODS: Twenty-seven healthy and 30 total knee replacement subjects (GMK Sphere, GMK PS, GMK UC) were assessed during multiple complete gait cycles of level walking, downhill walking, and stair descent using dynamic videofluoroscopy. Following 2D/3D registration, tibio-femoral rotations, condylar antero-posterior translations, and the location of the centre of rotation were compared. FINDINGS: The total knee replacement groups predominantly experienced reduced tibial internal/external rotation and altered medial and lateral condylar antero-posterior translations compared to natural knees. An average medial centre of rotation was found for the natural and GMK sphere groups in all three activities, whereas the GMK PS and UC groups experienced a more central to lateral centre of rotation. INTERPRETATION: Each total knee replacement design exhibited characteristic motion patterns, with the GMK Sphere most closely replicating the medial centre of rotation found for natural knees. Despite substantial similarities between the subject groups, none of the implant geometries was able to replicate all aspects of natural tibio-femoral kinematics, indicating that different implant geometries might best address individual functional needs.